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ABSTRACT 


The decision problem of the allocation of available weight and 
space in an air-dropped anti-submarine torpedo to fuel^ explosive^ and 
to attack speed capability above some preset minimum is studied for the 
case in which alternatives under consideration do not differ with respect 
to their effect on weapon reliability and to their effect on enemy 
countermeasures capability* Analytical models are developed which 
relate (l) probability of acquisition to endurance for the case of 
circling search at constant depths (2) probability of hit given acqui¬ 
sition to endurance and attack speedy, ( 3 ) probability of hit on a 
close-in attack run to attack speedy and (4) probability of kill damage 
given a hit to amount of explosive used* 
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CHAPTER I 


INTRODUCTION 

The Problem 

An air-dropped anti-submarine torpedo may be characterized by a 
long list of parameterso Specification of parameter values in a 
specific weapon is made such that the weapon may achieve a maximum 
of effectiveness in an actual tactical situation in wartimeo The 
state of the art in weapon design will determine the values of many 
of the weapon parameters* For example^ the transducer parameters 
may be prespecified by the status of some continuing torpedo sonar 
development program. Cost considerations may also be of great 
importance in the evaluation of alternative proposals for the 
design parameter values. For example, a choice between two propul^ 
sion systems may hinge on a difference in development costs necessary 
to produce a reliable weapon. However, for torpedoes and more 
especially for air-dropped torpedoes, the restrictions of weight 
and of physical dimensions (or volume) resulting from the varied 
mix of capabilities of the prospective weapon delivery systems are 
also very significant in the weapon design. 

Consider the restraints of weight and volume. One of the design 
parameters which has a great effect on weapon weight and volume is 
the maximum operating depth. An increase of a few hundred feet in 
the operating depth requirement will require a more substantial 
structure. This in turn may result in a decreased weight and volume 
to be made available for other weapon components. However, intelli^ 
gence reports or estimates of the depth capabilities of prospective 





I 







I 






submarine targets in the time period for which the weapon is to be 
designed may prompt the decision maker to fix the design operating 
depth as an input to the design problemo Similarly the advance sped 
fication of a minimum attack speed capability may be prompted by 
intelligence estimates and will certainly act to place an upper limit 
on weight and volume available for other weapon components« Suppose 
that a study of specifications from higher authority^ engineering 
capabilities to produce the various weapon components^ cost considera 
tions or some appropriate optimization procedures have resulted in 
specification of values for all weapon design parameters except 

1, Amount of fuel 

2, Amount of explosive 

Maximum speed (above some preset minimum value)o 
Further, suppose the restraints on weight and volume have been deter^ 
mined such that if the amount of fuel and explosive and the maximum 
speed capability above some specified minimum arc considered as a 
weapon payload, then the available payload weight and volume are 
provided to the decision maker* The decision problem of allocation 
of available payload weight and volume to amount of fuel^ amount of 
explosive, and to maximum speed capability above some specified 
minimum is the problem considered in this paper* 
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The Approach 

An initial approach to the allocation problem is to set certain 
minimum limits on the parameter valueso The minimum value of maximum 
(attack) speed has been specified« There must be sufficient fuel to 
provide some minimum probability of target acquisitiono (in the case 
of continuous circular search at a constant preset depth as is con¬ 
sidered in this paper a minimum probability of acquisition might 
result from a dive to search depth and completion of one complete 
search turn^) Then there must be sufficient fuel for at least one 
attack under some assumed conditions of acquisition range and target 
course and speed. Finally there must be at least enough explosive 
to result in kill damage if a hit occurs at the most favorable depth 
and point of contact on the target surface. Having determined these 
allocation minimums there remains the problem of allocating the 
balance of available payload weight and volume. 

A second step may determine a maximum value of attack speed. 

This maximum may be determined from a study including the probable 
operating depths the characteristics of the guidance system^ and the 
generated self-noise level and its effect on the ability of the sonar 
to maintain contact with the target on the attack run. 

What is the next step? How does the decision maker decide between 
enough additional fuel for one reattack and enough additional explosive 
to produce kill damage over a given additional element of target 
surface? What measure of effectiveness is used? 

This paper proposes a procedure which may provide the decision 
maker with a quantitative analysis or comparison of the tradeoffs 
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involved in alternatives of allocation among fuel^ explosive^ and 
additional attack speed capabilityo The tradeoffs are measured in 
terms of probability of kill as a measure of effectiveness<> Probabil-> 
ity of kill is dependent upon weapon reliability and upon the 
capability of the enemy to employ countermeasures (in addition to 
maneuvering)* It is assumed in this paper that weapon reliability 
and enemy capability to employ countermeasures (in addition to maneu^ 
vering) are both independent of the amount of fuel and explosive and 
of the weapon attack speed within the range of values between the 
preset minimum and some determinable maximum* This enables removal 
of reliability and enemy countermeasures capability (less maneuvering) 
from the problem and probability of kill may then be factored into 
three independent probabilities; those of probability of acquisition^ 
probability of hit given acquisition^ and probability of kill given 
hit* An optimization is made of probability of kill as a function 
of the fuel^ explosive, and additional attack speed capability with 
weight and volume restraints* This does not provide a number for 
use as an operational value of kill probability* However the 
dynamic programming methods used provide a solution which specifies 
the optimal allocation of weight and volume* 

It is emphasized that a quantitative comparison of the tradeoffs 
involved in the allocation process may be obtained only by some 
common analysis of the three study variables as they act upon the 
events of acquisition, hit, and kill* 
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The dynamic programming optimization model developed generates 
an input requirement for (l) probability of acquisition as a function 
of available fuel^ (2) probability of hit as a function of available 
fuel and attack speed and (3) probability of kill damage as a function 
of explosive weight. These requirements may be satisfied in a number 
of ways including 

1. analytical derivations in terms of obtainable or estimable 
inputs 

2« simulation (war gaming) 

3. extrapolation of empirical data from operational or 
development tests 

4. a combination of the above 

The model user may choose any of these methods so long as the 
resulting inputs are compatible with the allocation model. However^ 
any of the methods used will require an assumption about the tactical 
situation with respect to which the optimization is to be made. 
Consider the actual tactical situation in which the weapon will be 
used. Assume that methods exist to determine the optimum values 
of all weapon parameters, given that the weapon will be used in a 
completely specified tactical situation. The optimum parameter 
values determined will be different for each different tactical 
situation considered. The space of different tactical situations 
is infinite. However a single set of parameter values must be spe^ 
cified for the weapon, ultimately to be produced in quantity. By 
what criteria is the specific tactical situation to be selected for 
use in the optimization of the weapon parameters^ The probability 
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that the weapon will actually be used in the tactical situation selected 
is zero. The following two criteria will be used in this paper in the 
analytical derivation of the three required probabilities« 

Criterion 1: Select an unfavorable tactical situation which is not 
extreme. This may be thought of as a modified *maximin* criterion, 

A maximin criterion implies optimization of the weapon parameters to 
maximize effectiveness for that tactical situation which offers 
minimum probability of weapon success. The extreme situation is 
modified in this paper because the low probability of kill resulting 
from assumption of an unfavorable situation in the extreme would 
undoubtedly distort the relationships between the study variables which 
would apply in a greater percentage of cases. In general a weapon 
which is optimized for an unfavorable situation will be able to cope 
with a more favorable situation for which its parameter values were 
not optimized. 

Criterion 2: Select a tactical situation which is sufficiently 
simplified as to be within the scope of this paper. Progress in 
work of this kind is more often made by starting with something 
relatively simple and adding sophistication after experience is 
gained with the model. 
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Relation of this work to previous work 

The allocation problem has been studied before by R« Co Brumfield 
in references 1 and 2 in which algebraic expressions were derived 
relating total weight and volume of a weapon of arbitrary speed and 
range to performance characteristics anH component characteristicso 
This enabled determination of optimum running speed and range for 
stern chase vehicles under conditions of minimum propulsion system 
weighty minimum total weight of propulsion system and associated 
shelly and minimum propulsion system volume« The approach in this 
paper is essentially different from that of references 1 and 2 
although the problem is much the same« The previous work optimized 
total weight and volume and attack speed for the event of a single 
attack run, given explosive weight, a length of run in the acquisition 
phase, and acquisition range. This paper optimizes an allocation of 
payload weight and volume to fuel, explosive, and attack speed for 
the events of acquisition, hit (attack and reattacks) and kill 
given the maximum allowable payload. 

A great deal of work has been done in the study of various 
phases of the torpedo problem. References 3 4 form the background 

on which the development in Chapter III is based. Reference k 
develops a continuous volume scan rate for a single circular search 
turn. In this paper, in the first three sections of Chapter III, 
a volume scan for a single ping on the first search turn is developed. 
This is a modification of the continuous scan rate to the discrete 
case where the discrete interval is a single ping cycle. However, the 
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development in this paper is continued to obtain a probability of 
acquisition not only for a single circular search turn but for con«= 
tinuous turns until fuel exhaustion* This further development is 
made possible because a tactical situation is postulated in accordance 
with criteria one and two* The postulation is Chapter III limits 
the usefulness of the specifics of the development to the search 
pattern considered but is required in order to obtain the probability 
of acquisition input to the allocation model* The logic of the 
development may be modified as appropriate to provide for other 
postulated tactical situations in accordance with requirements to 
study a particular search pattern* The study of search turns subse¬ 
quent to the first and until fuel exhaustion is necessary to provide 
probability of acquisition as a function of fuel weight throughout 
the range of the argument * 

The development in Chapter IV of probability of a hit in any given 
pass at the target has no known precedent in the literature and is 
presented as an alternative way of looking at the hit probability 
problem which may be of interest* 

It is emphasized that Chapters III^, IVp and V are samples of 
possible methods of attacking the problem of finding inputs to the 
allocation model of Chapter I which are specifically designed to 
permit examination of the tradeoffs between alternative allocations 
of payload weight and volume among the three study variables* 
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Structure of the paper 

Chapter I contains the development of a mathematical model of 
the torpedo problem in a form such that the dynamic programming 
techniques of Chapter II of reference 5 applied to solve the 

allocation decision problem* Chapters III^ IV^ and V contain the 
development of the three required inputs to the allocation model of 
Chapter IIs probability of acquisition^ probability of hit given 
acquisition and probability of kill given hit* Appendix A is a 
description of the computer flow chart notation to be used in the 
flow charts of the other appendices* Appendix B is the computer 
flow chart of the dynamic programming allocation model of Chapter I* 
Appendices C 1 through C 8 contain the mathematical details of the 
developments in sections I through VIII of Chapter III plus the 
corresponding computer flow charts* This mathematical supplement 
to Chapter in is included in appendices to enable emphasis of the 
logic flow in Chapter III with a minimum of mathematical detail* 
Appendix D contains the computer flow chart for Chapter IV* 
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Summary 

A mathematical model is provided which enables optimization of 
the decision problem of the allocation of available weight and space 
in an air-dropped anti-submarine torpedo among three study variabless 
amount of fuel^ amount of explosive and attack speed capability above 
some preset minimum« The model is restricted to the study of alterna¬ 
tive allocations which do not differ with respect to their effect on 
weapon reliability and to their effect on enemy countermeasures 
capability (other than maneuvering capabilityo) 

Mathematical models are developed which relate (l) probability 
of acquisition to endurance for the case of circling search at 
constant depths (2) probability of hit given acquisition to endurance 
and attack speed, ( 3 ) probability of hit on a close-in attack run to 
attack speed, and (4) probability of kill damage given a hit to amount 
of explosive used. These four models are developed as specific 
examples of inputs needed to be compatible with the requirements of the 
allocation model. Computer flow charts are provided for each model 
except the last to facilitate programming in any desired computer 
language and also as an alternative method of study of the development 
of the model logic. 
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CHAPTER II 


THE ALLOCATION MODEL 

Consider the design of an air-dropped anti-submarine torpedo such 
that the amounts of fuel and explosive to be included and the capabi¬ 
lity of attack speed in excess of some predetermined minimum value 
may be considered together as a payload for which there is available 
a predetermined maximum of weight and volumeo Consider the employ¬ 
ment of the torpedo from the time it splashes into the water until 
the time it explodes in the proximity of a submarine target as an 
ordered sequence of the following three events each of which must 
occur in succession for accomplishment of the torpedo missions 

1. Acquisition after a search run following some preset pattern 

2. Hit after an attack (and if necessary a series of reattacks) 

Kill damage resulting from explosion of the weapon in the 
proximity ot the target 

Divide the study variable, amount ot tuel, into two separate and 
independent study variables: 

1. Amount ot tuel to be used in tne event of acquisition (event 1) 

2. Amount of fuel to be used in the event of hit (event 2) 

Restrict the study variable, attack speed, to a range of values 
such that 

1. The lower limit of the range is some preset minimum value which 
is greater than the given value of weapon search speedo 

2. As the attack speed is allowed to vary over its range, the fuel 
expenditure rate at attack speed and at any depth under considera- 
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tion may also vary. However as the propulsion system is modified 
to allow for changes in the value of its maximum (attack) speed 
capability, the fuel expenditure rate at search speed and depth 
remains approximately constant. 

Make the following assumptions^ 

Assumption 1: Weapon reliability and target capability to employ 
countermeasures other than maneuvering are both independent of changes 
(within the ranges to be considered) in the values assigned to the 
four study variables: 

1. Amount of fuel for event 1 

2. Amount of fuel for event 2 

3. Amount of explosive 

4. Weapon attack speed 

Assumption 2: The weapon is 1005^ reliable. 

Assumption The enemy target submarine has no countermeasures 
capability except that of maneuvering. 

Note: Assumptions 3 ^nd 4 are made in order to justify the relation^ 
ship between the study variables and the measure of effectiveness to 
be introduced later. The fact that they are intuitively incorrect 
does not detract from the validity of the optimization in the model 
because of assumption 1. 

Assumption 4: The three events of acquisition, hit and kill are 
independent events in a probabilistic sense. That is, if 
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P ” Probability that event x occurs 
= P ptarget acquisition (event 1)J 
Pjj = P ^target is hit (event 2) given that event 1 has occurredJ 
Pg = P jkill damage is inflicted on the target by explosion of the 
weapon warhead (event 3) given that event 2 has occurred 


If Pj^ = P [target kill] 

then 

Define: 

W = available payload weight (lb) 

V = available payload volume (ft^) 

= fuel weight for acquisition^ event 1 (lb) 
Wg = fuel weight for hit, event 2 (lb) 


W2*= weight for modifications to the propulsion system to provide 
a maximum (attack) speed capability in excess of some preset 
minimum value (lb) 

= preset minimum value of attack speed (yd/sec) 
w^ = explosive weight (lb) 

o 

= fuel volume for acquisition, event 1 (ft^^) 
v^ = fuel volume for hit, event 2 (ft^) 

v^*^ volume for modifications to the propulsion system to provide 
a maximum (attack) speed capability in excess of (ff^) 

O 

v^ = explosive volume (ft^) 

Wg = payload weight for hit, event 2 (lb) 

■ . * I 

W2 = ”2 ”2 
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^2 ~ payload volume for hit, event 2 (ft^) 

V, = V' . V' ' 

A W = increment of weight of the any of the weight variables, all of 
which are varied in increments of equal size (lb) 
gi(Wi) = the relationship between the i th weight variable and 

I V 0 

the i th volume variable^ i - 1^ 2^ 2 ^ 3 

III 

8i(w^) = v^j i = 1, 2, 2 , 3; g^(w) = g^Cw) 

fl(wi) = natural logarithm (in) of expressed as a function of Wj 

w** ) = In P expressed as a function of w* and w* ’ 

= In P expressed as a function of w^ 
j j f 3 

3 

F = In P = f 

i = 1 ^ 


I i 

W > W = lower and upper limits respectively on the allowable range 
of fuel weight (lb) 

II I 

W , W = lower and upper limits respectively on the allowable range 
of weight for modifications to the propulsion system to 
provide a maximum speed capability in excess of (ih) 


4 k • t t I 

w w 

^ ^ = lower and upper limits respectively on the allowable 

range of explosive weight (lb) 

The concept ot the lower and upper limits on the ranges of the 
study variables is as follows. It will be obvious to the design 
engineer that very small values or either fuel or the explosive 
weights will result in a negligible probability of success. Thus 
there is no need to examine the vicinity of the extreme values of 
zero and W. However the values of the limits may be chosen as 
close to zero and W as is desired. In the case of W*^ the value 
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chos6n may well be zero as this implies that the considered minimum of 
attack speed is to be the preset minimum« In the case of W the 
value chosen will probably be well under W as the maximum of the range 
of values of attack speed is limited by considerations of resultant 
guidance problems, excessive self^noise levels, and so forth« The 
hit event, event 2, contains two study variables which must be con¬ 
sidered together. For events 1 and 3 only one study variable is con¬ 
sidered per event (or programming stage). Thus for each value of 
fuel weight (w^^) there is associated a single value of acquisition 
probability P*(w, ) and a single value of volume v, =s gj^(wj^). 

Similarly for values of explosive weight w^. 

However 

+ w^' 

Tj = Tj « Vj' = gj(wj) * (Wji ) 

This requires that for each value of weight (w^) there be an 
associated set of values of hit probability volume 

V . one value for each combination of values of w' and w'* which 

d d. 

sum to the considered value of w^. A suboptimization of w^ and w^ 
to provide a single optimal value of PjjCwg) for each w^ as an input 
to the allocation model would neglect the associated volume 
restraints. Therefore index the weight w^ by k(k « K) over 

its range 

min -^W' , W'*^ — Wg —min -^W* + W % 
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and define 


w ~ ^ value of weight to be distributed between w^ and 


w^'; k = I,-.-, K (lb) 

'^2jk “ value of w^ associated with j - 1;,°°”;, Jj^ (lb) 

w^jj^ = j th value of w^' associated with j = 

^Hjk ” j value of associated with w^j^l j - 1#°°°# 

v^.^ = j th value of v^ associated with J “ °°°> 


2jk J — - 2 

k = 1, K 


^2jk “ ^ value of v^' associated with w^j^j j = 1 


» ° ° °» 'Ik^ 


V = V* , + vV., = j the value of v_ associated with 
2jk 2jk 2jk 2 

W2k; j = i> \ (ft^) 


The available payload W and the upper and lower weight limits or 
the weight increment A w are so chosen that all the weights are 
divisible by Aw. Then there will be 


W» W' 
Aw ” Aw 


considered values of Wj^ denoted Wj^^i r = 1, 
there will be 


min { W‘ + W» •, W J- min W», W> 


Aw 


Aw 


considered values of w^ denoted w^j^j k = 1, 


= K 

K 
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and there will be 


Will w" ’ 


Aw 


Aw 


= S 


considered values of w^ denoted w. t s = 1 , °, S 

3 3s' ' ^ 

The functions gj^(wj^), ggCwg), gg'Cwg')^ g^Cw^)^ and 

£3(^173) are assumed to be non-decreasing. That is^ an increase in 
the argument will not decrease the value of the function. In the 
case of '*^2* ^ assumed that any increase in either of 

the arguments while maintaining the other fixed will not decrease 
the value of the function. 

The problem is to maximize F = In 
subject to 

+ W3 

Vi + + V3 <V 

Consider the new problem to maximize 
P* = F - x|^gj(w^) + g^(wj) + 


subject to 


Wj^ + Wg + W3 <W 


10 


gi(w^) — 0 ; 1 = 1, 2, 2 i, 3 

where X is a non-negative real number (Lagrange-multiplier), 


Define the recurrence relations] 


h^(w) = fj(w) - X gj^Cw)' if W s w^^? r = 1, 


R 


= - 00 


otherwise 
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C ^ h^(W=W2j^) 


where 


^2^“2Jk’ "yk^ ^ %Jk ■ ^ *^B^"2Jk’ "yk^ 

’'ajk = e2("2Jk> * (wjj^) 


3(W) = ^ jfjCwj^) . XgjCwj^) ♦ hj(W.W3^)j 


The dynamic programming methods of Chapter II of reference 5 
next applied to the recurrence relations h^;, h^, and h^ for some fixed 
X to obtain a solution for the optimum value of h^(w) s f’* and for the 
corresponding optimum values of ^ ^3^ 

Next the corresponding optimum values of v^^^ v^, and v^ are com«=» 
puted using the computer determined optimum values of ^ ^3‘ 

The sum (vj^ then compared with the available payload 

volume V. Probably ^ ^3 ^ 


18 








Then if 


V 

the value of X must be increased and the dynamic program repeated 
until the minimum value of X is found which results in 

Vi + <V 

If v^ + Vg + v^ <V 

the value of X is decreased until the minimum value of X is found 
which results in 

Vi + Vg + <V 

The optimum value of h2(w), Wg* and determined for 

such a X specify the solution values of the three study variables® 
The Lagrange-multiplier, X > may be thought of as a cost 
coefficient® If: X = o 

Vi + Vg + <V 


then the volume restraint may be thought of as having cost nothing® 
That is, the measure of effectiveness has been optimized solely with 
respect to the weight restraint and the result has met the volume 
requirements. On the other hand, if X - o 

Vi + Vj + V3> V 

then some positive value must be assigned to X which will result 
in meeting the problem volume restraint Vj^ + v^ + v^^V. 

Ultimately the subsequent process of optimization with some 
positive value of X will result in smaller optimum values of one 
or more of the weight variables with a correspondingly reduced value 
of the measure of effectiveness. 
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This reduction in the solution value of the measure of effective- 


ness may be considered as the cost of having a volume restraint» Of 
course if weights and volumes are interchanged in the model the 
resulting X will be a measure of the cost of the weight restraint 
If the solution value of probability of kill is desired then 



" exp 
where 


, v^ and v^ are solution values* 


3 


The computer flow chart to accomplish this dynamic programming 
procedure is contained in Appendix B* 

It may be of interest to allow the values of the available 
payload weight and volume (W and V) to vary in the vicinity of their 
preset values* Then the sensitivity of the effectiveness measure^ 
to such small variations may be examined in order to determine 
the possible value of modifications to the weapon delivery systems to 
allow for a different (more effective) load capability* 

Input Requirements for the allocation model 
W Available payload weight (lb) 

o 

V = Available payload volume (ft*^) 


W ^ W = lower and upper limits respectively on the allowable range 
of fuel weight (lb) 

W , W = lower and upper limits respectively on the allowable 


range of weight for modifications to the propulsion 
system to provide a maximum (attack) speed capability in 
excess of a minimum preset value (lb) 
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W , W = lower and upper limits respectively on the allowable 

range of explosive weight (lb) 

11 

gl(wi); i =: 1, 2 , 3 = the relationships between the i th weight 

variable and the i th volume variable| 

I I 

i = 1. 2 , 3 

P^(wi) = probability of target acquisition expressed as a function 

of the fuel weight for acquisition 

P„(w*, w*’ ) = probability of hitting the target given that target 
rl 2 d 

acquisition has occurred expressed as a function of the 
fuel weight for hit and of the weight for modifications 
to the propulsion system to provide a maximum (attack) 
speed capability in excess of 

P (w ) = probability of target kill given that a hit has occurred 
expressed as a function of the explosive weight 
= minimum preset value of attack speed (yd/sec) 

^w = size of the increments of weight of the weight variables^ 
all of which are varied in increments of equal size (lb) 

The payload weight and volume are determined by the capabilities 
of the weapon delivery systems« The limits of the weight variables 
are determined by a preliminary study which reflects limiting design 
considerations and reduces the possible ranges of the weight 
variables to reduce computer timeo The weight^volume relationships 
are determined by engineering considerations reflecting the state 
of the art. The minimum attack speed is determined by intelligence 
estimates or CNO requirements. The size of the weight increments 
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is determined by considerations of required accuracy of results and 
available computer time. The probabilities of acquisition^ hit and 
kill may be obtained in a number of ways including 

1. analytical derivations in terms of obtainable or estimable inputs 

2. simulation (war gaming) 

3 . extrapolation of empirical data from operational or development 
tests 

k. a combination of the above 

The remaining chapters provide samples of possible analytical methods 
of providing the three required probability inputs which are specific 
cally designed to be compatible with the allocation model requirements 
It is possible to extend the model to provide for cost restraints 
The procedure would require a second Lagrange^multiplier in the 
function to be maximized (f* ), Additional inputs required would be 
the relationsuips between the costs and the corresponding weight 
variables. It is considered, however, that cost restraints are not 
significant to the problem considered in this paper. Variations in 
the amount of fuel and amount of explosive have little effect on the 
costs involved. Even changes in attack speed within the limits 
appropriate for consideration will probably have negligible cost 
effect when compared to the total cost of developing, producing and 
operating a weapon in the course of its life cycle. 
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CHAPTER III 


PROBABILITY OF ACQUISITION 

Consider a tactical situation in which a weapon is delivered by 
some air delivery weapons system to a splash point on the surface of 
the water in the vicinity of a submarine target and make the following 
assumptions. 

Assumption 1; The splash point defines the center of the top of a 
right circular cylinder called the target volume cylinder. The 
probability is quite small and is assumed zero that the target is out= 
side of the target volume cylinder. The target is placed at random 
within the target volume cylinder. 

Assumption 2: The weapon dives to its assigned search depth on a 
straight line track and commences a right circular search turn at 
constant depth, speed, and turn rate immediately upon reaching 
search depth. The weapon searches in the active acoustic mode 
commencing its first ping upon commencement of the search turn. 

The weapon continues its turn until target acquisition or fuel ex^ 
haustion. 

Assumption 3! The weapon does not acquire a false target. 

Assumption i<-: The search turn of the weapon is approximated by 

rotating the weapon about a vertical axis through the point at which 
the weapon commences search. 

Assumption 5; While the weapon is diving to search depth, the radius 

of the target volume cylinder increases at a rate equal to the target 

cruise speed. 
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Assumption 6: The target is alerted on commencement of the weapon search 
turn by the first 'ping' of the weapon. Thereafter the target opens at 
constant depth and speed radially from the projection of weapon position 
on the horizontal plane through the target. 

Assumption ?: During search the radius of the target volume cylinder 
increases at a constant rate equal to the target speed. 

Assumption 8; The 'definite range' law as developed in Chapter II of 
reference 6 applies to acquisition. Thus on a given radial during a 
ping cycle under specified conditions there is a range beyond which 
acquisition probability is zero and within which acquisition probability 
is one. Regions closer than the 'definite range' limits are considered 

scanned« 

Assumption 9: During search the weapoi makes a series of independent 
glimpses for the target, one for each ping cycle. 

Assumption 10; Target cylinder volume is considered constant during 
each ping cycle at the value valid at the initiation of the ping. 
Assumption 11; The water medium is considered to be isotropic. 

That is, all environmental parameters that affect acoustic performance 
have a constant value throughout the considered volume. 

Assumption 13; Consider planes containing the transducer and angles 
measured at the transducer. Consider the vertical plane containing the 
weapon longitudinal axis as the vertical reference plane. As the 
weapon sweeps through an angle <P in the horizontal plane it sweeps 
through a lesser angle <#>» in a plane perpendicular to the vertical 
reference plane but inclined to the horizontal. Also if the inclination 
of the plane inclined to the horizontal as measured in the vertical 
reference plane is « then the inclination of the plane inclined to the 
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horizontal as measured in a vertical plane other than the reference 

plane is a lesser angle . In this paper it is assumed that 

01 = 0 

oct = <x 

In practical situations likely to be considered, this assumption 
results in much smaller errors than those created by treating the 
range problem in discrete steps as is done in this chapter. 

Define; 

P sp (l) = Pr acquisition in 1 ping cyclesJ 
A A L 

where ^ “ probability of event X happening 

g. = P [acquisition during the i th ping cycle] j i = 1, I 

= weapon fuel expenditure rate at search depth and speed (Ib/sec) 

PRR = weapon active sonar pulse repetition rate for search 
s 

(pings/sec) 

w = fuel used by weapon to reach search depth, (lb) 

OSD 

Wj^ = fuel available for acquisition, (lb) 

Then, as developed in Chapter 2 of reference 6, 

(1-gi) 


and since 


i=l 


" ''OSD^ 


w 


It 


= I 


and |xij = greatest integer —X and where I 

is the number of complete ping cycles available for the event of acquisi= 
tion. it may be said that P. is a function of w^^ and a table may be 
computed of P^ vs where w^ is allowed to vary between W and W . 
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Of course if I is less than 1 then is defined to be zeroc The 
probability of acquisition in a specified glimpse {g^) may be con¬ 
sidered as a coverage factor <» The volume containing the target with 
positive probability is a region to be covered© A portion of the 
region to be covered is scanned (covered) by the weapon in accordance 
with the definite range law* Because the target has been randomly 
placed in the target volume cylinder (region of target uncertainty 
to be covered) then 

volume scanned on i th ping cycle 
^i " target volume during i the ping cycle 

Thus the problem is reduced to finding what portion of the volume 
of target uncertainty is scanned by the weapon on any ping* This will 
be different for every ping. 

The following development for is divided into eight consecutive 
steps. The development is amplified as appropriate and complete flow 
charts for a computer program are provided in appendices C 1 through C 8® 
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1. Available Echo Excess 

Consider a vertical plane containing the longitudinal axis of the 
weapon at the time of ping transmissiono Considerj, as in Figure Ij, a 


second plane perpendicular to the first and containing the weapon 
transducer such that it makes an angle j with the weapon longitudinal 
axis, A set of such planes is considered such that the angle j is 
associated with the j th plane. In the vertical plane angles are 
measured positive counterclockwise from the reference axis. 

Design or operational considerations will make it necessary to 
prohibit reception of signals from off the transducer axis in the 
vertical plane in excess of some limting angle J, not greater than 
90 degrees. Allowing one degree between each of the considered planes 

the index j will assume the range of values =J + 1, -‘“i) 

J. Thus 2J+1 planes are to be considered, one containing the trans= 

ducer axis and J planes on either sidep 


Consider increments of target range of 100 yards in the plane j. 

For an assumed target positioned in the j th plane at a range of 100k 
yards on the transducer azimuth axis projection in the plane j (point P 
in Figure l), there exists an echo excess by which the signal level 




of the echo reflected from the target exceeds the background 


masking level at the weapon transducer (m(BG)j^). That iss 


>09 If 
S> ^ 


Assuming a definite range law of acquisition, a negative echo excess 
yields zero probability of acquisition and a positive echo excess yields 
a probability of acquisition of one. The maximum range of acquisition on 
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Figure 1; 


The plane j. The transducer is at the origin pointing 
X direction. Line ox is the transducer axis extension. Lines 
ox, oP, and oy are in the vertical plane. Lines oP and oz are 
in the plane j. The polar coordinates of the point P in 
vertical plane using the x axis as a reference are 
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the azimuth axis in the j th plane is thus the range of lOOK^ yards at 

which B equals zero. At ranges less than lOOK. yards there is an 
jk . J 

echo excess of B., above that required for acquisition on the azimuth 
jk 

axis. This excess is available for use off the aximuth axis. 

The signal level varies with assumed target range and with the 
elevation angle of the target above or below the tranducer axis in 
accordance with 

SLjj^ = WS + TS - 2 + 10 log [^b^(j) 

where: 

WS = on-axis source level (db at 1 yard from transducer) 

TS = target strength of target (db) 

H = one way transmission loss at a range of 100k yards (db) 

Ic 

b^(j) = normalized vertical pattern function for transmitting 
pattern (i.e., the ratio of intensity at the angle j 
off vertical axis but on the azimuth axis to intensity 
on-axis in both the vertical and aximuth planes. 
b^(j) = normalized vertical pattern function for receiving pattern. 
There are several different kinds of masking noises and it is their 
total effect which is considered as the interference background masking 
level. Included are self-noise, target noise, ambient noise, and 

volume and boundary reverberation. 

There is also an enabling level designed into the weapon electrical 
circuitry to protect the acoustic system from excessive interference 
and to avoid interpretation of the interference as a signal. As used 
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here, the weapon enabling level (WEI^) includes the weapon recognition 
differential. It is a function of range (or time) and decreases with 
increasing range. The background masking level (m(BG)j^) is the larger 
of the two: interference background masking level and the weapon 
enabling level. 

The output of this step is a table of the available echo excess 

(V’ 
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Scanned Volume Radials; General Case 
In each of the planes j there is now defined for each value of 
range (100k) an available echo excess for acquisition on the 

azimuth axis. This available excess may be used to find limits for 
off-axis acquisition. Define an angled to be the angle through 
which the weapon turns while a ping travels to a range of 100k and 


returns, 
Then 




200 


where 

s weapon search turn rate (deg/sec) 
c = velocity of sound in the medium (ft/sec) 

For .cho ranging In tha plana J to a range o£ 100k, an energy pulse 
transmitted at angle p will return at the angle? -0^. 

The echo ekcass (p) reguired to acquire the target at range lOOk 

in plane j at angle p off the aximuth axis is 

A^j^(p) = - 10 log b^ (p) - 10 logbj^CP j^) 
where b^ ( p) = normalized aximuth pattern function for transmitting 

pattern 

( P) = normalized aximuth pattern function for receiving 
pattern 

If the required echo excess ( p) exceeds the available echo excess 

then acquisition will not take place at range 100k at an azimuth ; 
removed p degrees from the reference axis of plane j. However if 
A^k( P) does not exceed acquisition takes place. 
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jn each plane j and at each target range 100k there are two limiting 

angles between which acquisition will occur with probability one and 

outside of which there will be no acquisition. Looking the trans=. 

ducer toward the potential target define the set of limiting angles 

6 with m = 1 for the limiting angle to port and m = 2 for the 
^jmk 

limiting angle to starboard. As shown in Figure 2^ the reference 
axis for polar coordinates in the j th plane is the projection of 
the transducer axis at the time of ping transmission in that plane. 
Positive angles are measured clockwise from the reference axis in all 
planes j and in any horizontal plane. The pattern may not be centered 
on the axis and thusp^^j^ may be positive or negative. Design or 
operational considerations will make it necessary to prohibit reception 
of signals from off the transducer axis in the azimuth plane in excess 
of some limiting angle * in the direction opposite to the search turn 
and er ' in the direction of the search turn. This restricts the 
angles P^^j^ as follows 

^jlk ‘ "^k 

Pj2k " k I “ * 

Recall that in plane j the maximum range for a non-negative 
was computed to be lOOK^ yards. For a straight running weapon this 
would be the maximum acquisition range in plane j. Actually an 
acquisition range of IOOk'^ yards will not be realized in plane j in 
this problem because the weapon is turning at some rate. The energy 
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Figure 2t Acquisition limiting angles in the plane j. Line OP is the reference axis 
The X axis is the transducer axis projection. In this c&s &^zero 




















from a given source direction will never return as an echo in the same 
direction. Thus there will always be some penalty for off-axis trans¬ 
mission or reception. 

Unless the turn rate is extremely high, the insonified volume of 
adjacent weapon pings will overlap. If a given ping scans part of the 
volume insonified by the previous ping, the overlap cannot be credited 
to the volume scanned by the later ping because the overlap volume was 
removed from the volume of target uncertainty by the previous ping. 

The concept of increased probability of detection because of overlap 
is not valid in this case because of the use of the definite range 
law. Thus all overlap volume should be credited to the ping which 
first scanned it. This would give all pings subsequent to the first 
ping during the first search turn an equal scanned volume less than 
that of the first ping. In this paper the situation is approximated 
by dividing the volume of overlap between the pings concerned to give 
all first turn pings an equal volume of insonification. This approxima¬ 
tion decreases P (1) but the error in P (l) decreases with increasing 

A A 

I until the end of the first search turn when it becomes negligible. 

The plane of maximum overlap is plane 0. In this plane, for small 
k, pQ 2 j^ ^or ping i will be larger than ^ P pi^g (i+l)« 

See Figure 3* Here h is the angle between the axes of pings i and 
(i+l) or the angle through which the weapon turns between pings. As k 
becomes larger the differences between P^^j^ for ping i ^ P 
ping (i+1) become smaller. Define p* to be the angle at which the 
two are equal. Thus the right-most angle for acquisition in any plane 
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Figure 3; Overlap of insonified volume between adjacent pings. 
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is not greater than p' and the leftmost angle is not less than p = b 
Some of the extreme planes j may not be useful for acquisition at all 
By including for consideration only those planes in which acquisition 
may take place the index j takes on the new range from j' below the 
axis to j'' above the axis. 

The range or radial distance associated with the angle is 

100k yards. See Figure 4. There is a range maximum called the ping 
interval range defined by the time between pings (l/PRRg) and the 
velocity of sound in the medium (c). Further, for planes inclined 


above the horizontal there is another range maximum defined by the 
range at which the radial at angle breaks the surface of the water 

Similarly for planes inclined below the horizontal there is a range 


maximum defined by the range at which the radial at angle p^^j^ breaks 
the maximum target depth (d^^). Volume scanned below the maximum 
target depth does not contribute to acquisition probability. Thus the 
acquisition range at any angle p.^j^ is the minimum of the ranges just 
discussed to be denoted 

Since the index k is no longer of value to indicate range it is 
useful now to relabel the radials such that the index k indicates the 

separation from the leftmost radial in the plane. See Figure 5- The 

leftmost radial is labeled (R.^, p^^). The general radial is 

(R Pi )• The rightmost radial is (R.j^ » PjK.^' planes close 

' jk' ^jk -• j ■’ J 

to the axis, probably 


PjKj = P ' PjK. ■ Pjl " ^ 

This may not be true of the extreme planes. 
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Figure k: Acquisition range at angle Line OR has length 

look. Line SS' denotes the surface. Length 01 is the 
ping interval range. Lined area is scanned. 
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Figure 5: Relabeled radials in plane j. 
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The output of this step is tables of 

V V’ ^ * ■’’’ ■■■’ ''“ 

plus the angles p' and b . 


K. 
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3* Scanned Volume: General Case 

The radials computed in step 2 define areas of insonification in 
each of a set of planes j. The area in each plane may be computed as 
the sum of a set of triangles as shovm in Figure 6. 

Aj = (1/2) gp Rjk ^j(k+l) 


The irregular shape of the areas in the planes j makes it difficult 
to work with them in their present form. Define a set of isosceles 
triangles (one for each plane j) such that the triangle areas are 
equal to the scanned areas A. in planes j and the common apex angle 


is the angle turned by the weapon per ping (?)). 
altitude of each triangle is 

r IV2 


H 


i 


tan(6/2) 


See Figure 7. The 


The set of triangles so defined approximates the areas scanned and 
makes it possible to obtain the total scanned volume SV^j^ on ping zero 
on the first turn. Ping zero is a ficticious ping denoting the 
general case. The computation of (scanned volume on i th ping of 

first search turn) is identical to that for to be developed next. 

Define; p= weapon angle of attack (pitch) (deg). 

The plane j makes an angle of (j+p) degrees with the horizontal. 

The set of isosceles triangles is used to define a set of pyramids. 
A vertical cross-section of the set of pyramids is shown in Figure 8 
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Figure 6: Area of insonification in plane j. 
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Figure 7 : 


Equivalence of area and isosceles triangle 
OHH* of altitude and apex angle b . 
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ping zero, x is the horizontal axis, y is the 


vertical axis, and OQ is the transducer axis 
extension. 


1^3 







and a single pyramid is shown in Figure 9- The volume of a single 
pyramid is 

Vj = 1/3 ^area of basej j^altitude^ 

As shown in Figure 9> 2X is the average value of one of the 
base dimensions and Y is the value of the other dimension. 

Thus the base area is 2XY. 


X = ||(Hj + tan(6/2)]/2 

Since sin0 = h/tlj = sin 1 )/l 


the altitude is 

h = (Hj sin 10 )A 

and 

Vj = (l/3)(2XY)h = sin 1® tan(&/2) 


The volume scanned is 


J"-l 


SV«, = ^ 


'01 


j=J' 


The output of this step is tables of j = J*# j'* 

and the volume SV^j^. 
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4. Special Volume Elements: General Case 

For the case of search at constant depth in a continuous circle 
as is being considered in this paper, the situation changes drastically 
upon completion of the first search turn circle <> Subsequent search 
turns will effectively scan much less volume because the weapon retraces 
its path of the previous turn. Define the maximum available number of 
search turns (u) to be 


U = 



"oSD 


) PRR„/^ 
s 



+ 1 


where M is the number of ping cycles in one search turn 


M= ||360/6|| +1 


l-j- 


= greatest integer x 


x|| = greatest integer <; x 


W = maximum allowable weapon fuel (lb) 

Wj^^ = fuel expenditure rate at search depth and speed (ib/sec) 
Wqsd = fuel used in dive to search depth (lb) 


°WS ''ot 


w, 


OSD sin<;£.jj5 


1/3 = weapon search depth (yd) 

w ^ 

= average weapon fuel expenditure rate in dive to search 
depth (Ib/sec) 

= weapon speed in dive to search depth (yd/sec) 

0 = angle of dive to search depth (deg) 

i/D 
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Consider a side view or vertical cross-section of the scanned volume 
of a single ping as in Figure lOo The vertically lined area is that 
scanned by the weapon on the i th ping of the first search turn,, 

Recall assumption 6 that the target opens radially from the weapon 
position and assumption 12 about weapon orientations on different 
search turns. Define a distance of TES yards as the distance traveled 
by the target during the time it takes the weapon to complete one 
search turn. The weapon scanned area on the i th ping of the second 
search turn is equivalent to that of the areas lined sloping downward 
to the right on Figure 10. The weapon actually scans the same 
geographic volume on each search turn but since the target is assumed 
to open radially the situation is equivalent to translating the 
area scanned on the first turn a distance TES to the left for each 
subsequent turn. If the target is in the area lined sloping downward 
to the left on the i th ping of the first search turn^ he will escape 
prior to the i th ping of the second search turn. Similarly for sub« 
sequent search turns. A target above the critical depth CDA or below 
the critical depth CDB may escape but a target at any intermediate depth 
will be acquired eventually if fuel supply permits because of the assump- 
tion of an opening target. The search turns after the first are desirable 
because of the possibility of acquiring a target * trapped® inside the 
insonified volume of the first search turn. Treatment of scanned 
volumes on the second and subsequent turns will differ depending on 
whether the volume is above or below the depths CDA and CDB. Although 
the areas scanned above depth CDA on search turns 1 and 2 are equal in 
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Figure 10s Vertical cross=section of a single ping. Lines xx* have length TES 
Lines yy® and lines xx* have length L . 
















































































Figure 10^ the volumes represented by these areas are not equals A top 
view of ping i is a pie shaped slice which makes the volume associated 
with the first search turn area larger than that of the corresponding 
second search turn area. Thus each volume element must be calculated 
independently and the volume scanned by the weapon on ping i of search 
turn u is the sum of several independently calculated volumes. 

In this step the volume elements in the general case (i-O) are computed. 
These are modified as necessary for each ping orientation in a later 
step. Although Figures 11 and 12 and the development to follow involve 
the volumes above the search depths the development for volumes below the 
search depth is identical. 

Consider the vertical cross-section of Figure 10. The cross- 

section has been split into several regions or volume layers identified 

on the left margin by their bottom boundary labeled such that the 

bottom boundary of a level (or layer) at or above depth CDA has a 

length of L, inside the first turn scanned region. The planes j are 
a 

considered using the plane (j = TJ) of maximum |sin (j + p) j as a 
reference plane. In Figure 10 below search depth 
I sin (TJ + ji) I = (TJ+l+^i) | . 


In this case where more than one of the planes j has the same maximum 
I sin (j + p)|, the reference plane is chosen to be that with the 
least cos (j + p). Above search depth, each plane j such that 
H. sin (i + p)<D,,- - CDA and i^TL defines a level L,. Also each 
plane j such that sin (j+ CDA defines a level 
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Referring to Figure 11^ the volume elements LADV^^;, MADV^^^ and 

RADV^^ (d = 1, D; u = 2> U) are significant* The subscript d 

indicates the level and the subscript u indicates the search turn 

number. The level L- is the final level above the critical depth CDA 

and = TES. The corresponding volumes below the search depth are 

LBDV^^, MBDV^^, and RBDV^^ (d = 1> ""**> ^ 2^ U )* The scanned 

volume in level above the search depth on the u th search turn is 

ASV , = LADV, + MADV^ + RADV^ 

du du du du« 


One or more of these volume increments may be zero as is LADV^^ 

in Figure 11. The corresponding scanned volume below the search depth 

is BSV , . 
du 

Referring to Figure 12^ 

the volumes d = D+lj D*j u = 1^ are significant* 

The corresponding volumes below the search depth are denoted 

(d = D -f 1, d'; u = 1, U ) 

For levels L, below depth CDA and above search depth and for the u th 
a 

search turn, the trapped volume of target uncertainty remain^ 

ing after the ping cycle is complete* The volume scanned in level 
above the search depth on the u th search turn is thus 


ASV^ = AMV^/ , V - AMV , 
du d^u-l; du 

The corresponding scanned volume below the search depth is 
In this step all calculations are done for ping zero and the 
quantities computed for ping zero are saved to be modified later as 
appropriate for the i th ping. The total scanned volume for ping 
zero of the u th search turn is 
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D» D« 

sVou= 2 “I'au * S 

d=l d=l 

The output of this step is the number of pings (m) in a. search 
turn^ the maximum available number of search turns U<) several tables 
of parameters described in appendix Ck used in computing the special 
volumes, and tables of the following special volumes for ping zero; 


LADV_, , MADV_, , RADV, ; d = 1, 
du' du' du' ' 

D| 

u = 2, 

o • • ^ U 

“DVju- ““''du' ““''du! 

*• "> 

u = 2, 

• • ® ^ u 

= 1, •••, d’; u . 

1, •••, 

U 


BMVdu! = 

1 « • « 

U 


ASVj^J d . 1, •••, D'i a = 2, 

• • u 



®^''du' ■* = 2. 

• • u 
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5. First turn scanned volume: i th ping 

Throughout the search phase^ the weapon is displaced from the 
center of the target volume cylinder a horizontal distance The base 

radius of the target cylinder is some distance at the time the 
weapon splashes into the water. The radius then increases at a rate 
(Stc) equal to the target cruise speed to a distance at the time 
of the first weapon ping. Thereafter the radius increases at a rate 
(S equal to the alerted target escape speed and has a value R. 
upon commencement of the i th ping cycle. 

In the situation of Figure 13^ a top view of the target cylinder 
of increasing volume with the first few ping orientations and scanned 
volumes in the first search turn, the insonified volume of ping 1 
includes a large volume outside the target cylinder. Such volume does 
not contribute to acquisition probability because there is an assumed 
zero probability that the target is outside. Therefore the total 
insonified volume must be reduced to the useful scanned volume inside 
the target cylinder. In the situation depicted, the amount of dis« 
carded volume is reduced with each ping until ping 9 when all of the 
insonified volume is within the target volume cylinder. 

The projection of the radial plane j to the hori= 

zontal plane through the transducer is the line OP of Figure Ik- 
denoted 

where 

"jk' ®jk > 
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Figure 13: Top view of target cylinder of increasing 

volume with the insonified volumes of several 
sample pings of the first search turn depicted« 
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Figure 14; The projection of plane j in the horizontal at 

search depth. Line OO' is the plane j reference 


axis projection. Line IP is 
r'., . Line IP' is R, . Line OP' is r' 




Line OP is 
ijk“ 
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Use of the law of cosines enables determination of the distance 


I P denoted This distance is compared to the radius (R^) of the 

target cylinder on the i th ping. If the radius (R^) is smaller^ it 

defines the new triangle lOP* and the projection of the side OP^ back to the 

plane j defines a new acquisition range (r,.. ) to be associated with the 

i jk 

angle p,, . That is 
J I 

^ik 


R. 


ijk cos 




However if the distance CR., is the lesser of the two then 

jk 


R. = R,, 

ijk jk 


The remaining portion of the computation for scanned volume (SV^j^) 
on the i th ping of the first search turn is identical to that in 
step 3 for S^oi* 

The output of this step is the table of volumes i s 1^ ^ M“1 

and the distances R and R.: i = 1, M-1. 

o 1 
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6. Scanned volumes subsequent to first turns i th ping 

Just as the first turn scanned volume had to be corrected for the i th 
ping because of the effect of the target cylinder radius, so must the 
scanned volumes on subsequent turns• 

Define the intersection of the ping i bisector with the target 
cylinder radius at transmission of ping i to be a distance CH from the 
weapon• 

Referring to the vertical cross-section of the insonified volume 
in Figure I 5 , verticals through four of the possible intersections have 
been plotted and labeled CHI through CH4 o The effect of each of these 
four cases will be discussed as it applies to a level above and a level 
below the critical depth CDA (i,e: levels 2 and 4)o The treatment of 
the situation below the search depth is identical« 

In the case of CHI the first turn scanned volume is reducedo 
However the scanned volumes for all subsequent turns are the same as 
those computed in step four for the general ping* This is because 
CHI lies to the right of all areas representing equivalent volumes<> 

In the case of CH2 there is still no reduction for volumes in 
level 2. However there is a small reduction in the volume scanned on 
the second search turn in level The volume reduction is that 
volume associated with the area in level 4 to the right of CH3o 

It is computed as a reduction in (the volume missed in 

level above search depth on the first search turn)o 
Thus i equals but 

‘“\2)piu8 i ('‘=%2)plng 0- 
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as numberedo Escape volumes have primed numbers according to the search turn of escapeo 
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In the case of CH3 there is a reduction in the scanned volume on 
the second search turn in both levels and o In level the volume 
(ASV^^) . . is less than (ASV^^) . * by an amount equal to 

RADV 22 plus a portion of 

In level L|^ the reduction still amounts to a portion of 

In the case of CH4 there is no scanned volume in level above 
search depth on the second turn^ That is - 0 for ping i. 

Further there is a small reduction in level in the scanned volume 
for the third search turn* For level L|^ there remains some scanned 
volume on the second turn and there is also a reduction in scanned 
volume for the third turn. 

The output of this step is a table of scanned volumes 
SV^^; i = 1, •••, M-1; u = 2, U. 
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7. Target Volume: i th ping 

Consider that portion of the target cylinder volume computed at 
a base radius of within the radials bounding the ith pingo See 
Figure l6« The area contained between the boundary radials of the 
i th ping and within the cylinder of radius multiplied by the 
height of the cylinder is this volume* The probability that the 
target is in this i th volume slice of the target cylinder is the 
volume of the slice divided by the total cylinder volume* This 
probability i M-l) remains constant throughout the 

problem because of the opening target assumption number six* 

From Figure l6: 

PQ _ area of horizontally lined region 
^ TT 

The target volume (XV^; i = 1, M"l) associated with the i th 

ping on the first search turn lies between the same radial limits but 
within a target cylinder of radius instead of Rj^. That is 

TV^ = (area of vertically lined region) 

Figure 16 shows the angle (&' ) assigned to the M th ping to be less 
than6 . It is 

b ' = 360 - (M-l)b 

This smaller angle is required to complete the search turn so 
that the i th ping of the u th turn is oriented in the same direction 
as the i th ping in the (u + 1) th turn. Thus the M th ping is in 
part ficticious. 

The output of this step is thus the tables of TV^j i s 1^ ••<>^ M=1 
and of i = 1, '•*> M-l and the angle b* » 
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Figure l6: Top view of target cylinder of increasing volume . 

Lines CT and CT' are boundary radials of ping i. 

2 

Horizontally lined area is PQ^ times . 

Vertically lined area times cylinder height is 
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8 . Acquisition Probabilities; 1 pings 

The instantaneous probability of acquisition on the i th ping is 
simply the ratio of the volume effectively scanned to the volume of 
target uncertainty. 

Define; g^^ = P j^acquisition during the i th ping cycle of the u th 

search turn]] 

Then PQ^ / TV^ J i = 1, u = 1, •••, U 

Recall that PQ^ is the probability that the target is within the 
radial limits of ping i. 

The volume TV^ is the total volume of target uncertainty within 
the radial limits of ping i. The volume is the volume scanned 

on the first turn. The volumes scanned on subsequent turns are volumes 
measured inside the first turn scanned volume and under the conditions 
existing for the first turn. These volumes are then moved intact with 
the opening taget. Thus the applicability to all turns of the volume 
TV^ measured using the first turn cylinder radius. 

The instantaneous probability of detection on the M th ping is 
approximated by 

= 8(M-1)u * ei(u.i) 

The M th ping is alloted the same time interval as any other pingo 
Errors involved should be small. 
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There reniatns the requirement to relabel the to be 

g^; i = 1, to compute I as a function of as described in 

the beginning of the chapter, to compute P as a function of and 

ri. 

to tabulate the values of versus as is allowed to vary from 

I 

W to W . This final table is the input to the dynamic programming 
model. 
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9. Input Requirements 

A. Medium 

c = velocity of sound in water (yd/sec) 
a = attenuation coefficient (db/yd) 

10 log m^ = surface reverberation scattering coefficient (db) 
10 log m^ = volume reverberation scattering coefficient (db) 

B. Target 

TS = target strength (db) 

TN = radiated noise level of the target (db at 1 yd) 

= maximum target depth (ft) 

S ^ = target speed during weapon search (yd/sec) 

X D 

S„ = target cruise speed (prior to becoming alerted on 
XU 

the first ping cycle) (yd/sec) 

C. Weapon 

WS = transmitted on-axis source level of weapon (db at 1 yd) 
AT = ping length (sec) 

PRR = weapon active sonar pulse repetition rate tor 
s 

search (pings/sec) 

Af = receiving bandwidth (cycles/secj 

WELj^ ss weapon enabling level tor a range of lOuk yds (dbj 
SN = equivalent isotropic sel£=noise level at search 
speed, deptn, and turn rate (dbj 
RD = recognition differential for 50^ probability of 
acquisition (db) 
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=s normalized pattern function in the vertical plane 
for transmitting pattern (a numeric) 
bg = normalized pattern function in the vertical plane for 
receiving pattern (a numeric) 
s normalized pattern function in the azimuth plane for 
transmitting pattern (a numeric) 
b|^ = normalized pattern function in the azimuth plane for 

receiving pattern (a numeric) 

J = limiting angle for off-axis reception in the vertical 
plane (deg) 

or, or» = limiting angles for off=axis reception in the azimuth 
plane (deg) 

DR = directivity index for receiving (db) 

p = angle of attack (pitch) of weapon (deg) 

= search rate of turn (deg/sec) 

DS = angle of dive in diving to search depth (deg) 

D = weapon search depth (ft) 
w ^ 

= weapon speed in dive to search depth (yd/sec) 

S,,_ = weapon speed during search 
wS 


W = maximum allowable weapon fuel weight (lb) 

Wit = weapon fuel expenditure rate at search depth and 
speed (ib/sec) 

Wq^ = average weapon fuel expenditure rate in dive to 
search depth (ib/sec) 
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D. Delivery weapons system capability 

= radius of target volume cylinder at time of weapon splash 
Statistical averages and distributions are available for the medium 
inputs* Intelligence reports and estimates or CNO requirements may 
provide the target inputs except for the target speeds which may be 
selected in the higher ranges in accordance with criterion one* 

The weapon inputs are available from weapon designers at least in the 
form of educated estimates. Finally the input may be estimated 
from observations of the performance of available delivery weapons 
systems. 
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CHAPTER IV 


PROBABILITY OF HIT 

Consider a tactical situation in which a weapon acquires a submariiue 
target at some initial slant range and make the following assumptions» 
Assumption 1: The weapon makes an attack and a series of independent 
reattacks on the target culminating in either a hit or fuel exhaustion* 
Assumption 2; At the time of acquisition^, the target is opening at 
maximum speed* Target course and speed are constant throughout the 
attack and reattack* 

Assumption 3: On the attack and after reacquisition on each reattack 
the weapon maintains contact and closes at attack (maximum) speed* 
Assumption k: The target operates in depth within specified minimum 
and maximum limits such that the maximum depth limit is not greater 
than the weapon maximum operating depth* 

Assumption 52 The depth at which a hit may take place is a random 
variable denoted the hit depth which has an assumed probability 
density function* 

Assumption 6: The average weapon depth for the attack is the arithmetic 
average of the search depth and the average hit depth* The weapon and 
target depths for the final stage of the attack and throughout each 
reattack are equal at some depth within the limits of assumption four* 
Assumption 7 : The attack and each reattack is from astern the target* 
There is some close-in range at which the weapon sensors lose contact 
with the target* On losing contact with the target the weapon continues 
on a straight track at attack speed for a preset dead time* Upon com° 


68 








pletion of the dead time run, the weapon instantaneously slows to 
search speed and commences a left turn at reattack turn rate and 
constant depth. 

Assumptions: The weapon reattack turn is a circle with no advance or 
transfer. The weapon reacquires the target upon completion of one 
complete reattack turn and instantaneously increases to attack speed. 
Assumption At the time the close-in weapon sensors lose contact 
with the target, there exists an intersection of the weapon track 
extension with the vertical plane through the submarine target center 
which is perpendicular to any plane containing both the target and the 
weapon. Let the target projection in the vertical plane be approximated 
by a circle denoted the target circle centered on the target longitudinal 
axis with a diameter equal to the length of the beam dimension of the 
target. Let the circle center be the origin of a doubly infinite 
axis containing the intersection. Assume that the intersection co¬ 
ordinate on the axis is a normally distributed random variable with 
mean zero and standard deviation o • 

Assumption 10: A hit occurs if the weapon track at the closest point 
of approach to the target center passes within a distance of the center 
not greater than the radius of the target circle plus a small proximity 
correction factor. 

Define: 

Pr = Pr (N) = P l^hit in N attacks] 

= P ^hit on the n th attack] ; n = 1, N 

Pjj (0) = 0 
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Then 


N 

Ph = 1 - TT (1 - 

“ n=l ^ 

The derivation is analogous to that for acquisition probabilityo 
Define: 

t^ = time the weapon sensors lose contact with the target on an 
attack or reattack (zero reference time) 

RH = range from the weapon to the target center at time t^ (yd) 
t^ = time from t of closest point of approach of weapon and 

UirA O 

target (sec) 

MRH = distance traveled by weapon from time t^ to time (yd) 

S = weapon attack speed (yd/sec) 

WA 

= target maximum speed (yd/sec) 

BW = target beam width (yd) 

VT = weapon proximity correction factor (yd) 

r = coordinate of intersection of the weapon track extension 

with the vertical plane (a) through the submarine target center 
perpendicular to any plane containing both target and weapon at 
time t^ (yd) 

r* = coordinate of intersection of the weapon with the vertical plane 
(a* ) parallel to plane A containing both the submarine target 
center and weapon at time (yd) 

F^(x) = cumulative normal probability distribution with mean zero and 
standard deviation one evaluated at the point x 
Consider the attack situation depicted in Figure 17^ 

The xy plane is the vertical containing the weapon and target at time 
t^* Planes A and A* are parallel to the yz plane and contain the 
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Figure 17o Planes A and A* in the attack situation 






















bases of two right circular cones with a common apex at the origin« The 
line OR is the weapon track and its orientation with respect to the x 
axis determines the orientation of the axes LL and L® on which 
coordinates r and r* are measured. The common intersection of the cone 
of base radius r*, the plane A*, and the xy plane is a line of length 

2|r*| . By similar triangles 

I MRH 

^ RH ^ 

Since 


MRH = RH - Sj„) 


then 


®WA / 


Thus r* is a random variable with mean zero and standard deviation 0*^ 


^ ^WA / ^^WA “ 


Therefore 


h = 


P 1^2 |r‘ I BW + 2 Vt] 
= P Qr' I < (BW/2) + Vt] 


= 2 F 




(BW/2) + VT 
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(BW/2) + VT 

o / (Sy^ - S^J^) 


Define; 

HD = hit depth (ft) 

= probability density function of the random 
variable hit depth 
D^j^ = maximum target depth (ft) 

*^MIN = minimum target depth to be considered (ft) 

Define a set of Q equal depth subintervals in the depth interval 
^TmJ that the depth at the lower boundary of the q th 

°MIN ^ 

= P Qxit depth is in the q th depth subinterva^ j q=l;i °) Q 

Assume as in Figure 18; 


subinterval is 


Define 


- 2 (x “ / (D„ 


- Dj^)' 


Then; 



gq - 1 . 


q 


1 , 


Q 


Define; 

wl = fuel used in the n th attack or reattack (lb) 
2n ' 

= initial acquisition range (yd) 
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Figure l8: 



Assumed probability density function of the 
random variable hit depth. 
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= weapon fuel expenditure rate at attack speed and 
at the average depth for the attack, (ib/sec) 

HD = average hit depth (ft) 

Then 

“ ' 2/3 ^ 1/3 

and the average depth for the attack is 


WS 


+ HD 


“ 1/2 D„s + 1/3 Dtv + 1/6 D, 


TM 


MIN 


For the attack 


w 


W = 
21 


2t 


WA 


TM 


Define: 


w 


Stq 


w 


Atq 


'DT 


'CT 


% 


= weapon fuel expenditure rate at search speed and at 
the midpoint of the q th depth intervalo (ib/sec) 

= weapon fuel expenditure rate at attack speed and 
at the midpoint of the q th depth interval. (Ib/sec) 

= dead time from t^ at which the weapon commences the 
reattack turn (sec) 

= time from t^ at which the weapon completes the reattack 
turn and reacquires the target (sec) 

= reattack rate of turn (deg/sec) 


Then: 


t 


CT 


'DT 




t 


CPA 


RH 

^WA ■ ®TM 
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I 




I 




i 




The weapon will spend the following time at attack speed during each 
reattack; See Figure I 9 . 

~ ^TM~ ^^DT°^CPA^ \a 

A = DT - CPA ^ 

The weapon will spend the following time at search speed during 
each reattack: 

T - 

S= 0^ 

Therefore; 


Q 

”2n = 21 "At, * '^S "Stq ] % 

q=l 


n = 2, N 

Define; 

S = minimum preset value of attack speed (yd/sec) 
o 

w^ = available fuel weight for hit, event 2, (lb) 

Wg' = available weight for modifications to the propulsion 

system to provide a maximum speed capability in excess 
of S (lb) 

N is a function of w^ according to 


N 


N+1 


w. 


2 n 




w, 


2 n 


n=0 


n=0 


However, recall that w* is a function of attack speed which is it= 
' 2n 

self a function of w^'. 
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Figure 19: Geographic plot of weapon and target tracks 

separated to show time differences* The 
distancesbelow the target track are measured 
from the joint position at time t^^.• 

w xtA 
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Thus 


N(w^, w^') 

p„(N) = 1 - -FT- [i - h„(v,' ') 
n=l i- 


Recall Wg = Wg + w^' 

W, W = lower and upper limits respectively on the range of fuel 
weight (lb) 

» » f 

W ^ W = lower and upper limits respectively on the range of 
weight for modifications to the propulsion system to 
provide a maximum speed capability in excess of (lb) 
For each value in the range for w^ 
min -|^W + W , 

the function P is evaluated for each possible combination 
n 

of the variables w^ and w^* as they are allowed to vary such that 

they remain within their respective ranges 

I I 

W ^ w. ^ W 
— 2 

M II _l I 

W <W2^W 

and their sum is equal to w^o Let K be the number of values of 
Wg. This procedure provides sets of values 

2 jk^ 2ik^ Hjk^ ^ ^ k 

for the k th value of (k = 1, K) 

The entire set of sets of values 


^2V.* ''2jk" ''2jk^ ^Hjk' ^ 
k = 1, ” ” K 


o o o 

> 


J 


k 
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is the input to the dynamic programming allocation models 

Input Requirements 

1. Target 

BW = target beam width (yd) 

DtM “ maximum depth (ft) 

= target maximum speed (yd/sec) 

2. Weapon 

S = weapon search speed (yd/sec) 
w o 

VT = weapon proximity correction factor (yd) 

0 )^ = weapon reattack turn rate (deg/sec) 

I 

W, W = lower and upper limits respectively on the allowable 
range of fuel weight (lb) 

= minimum allowable weapon attack speed (yd/sec) 

II ^1 I 

W , W = lower and upper limits respectively on the range 
of weight for modifications to the propulsion 
system to provide a maximum speed capability in 
excess of S (lb) 

Gj = Relationship between increased weapon speed above and 
increased weapon weight 

= Relationship between fuel expenditure rate and speed and 
depth. Recall that this relationship may vary as a 
function of attack speed only in the range for speeds 
above search speed. This relationship is required to 
specify and (ib/sec) 
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3. Miscellaneous 


RH = range from the weapon to the target center at the instant 


(t^) during close-in attack when the weapon sensors lose 


contact with the target (yd) 

0 = standard deviation of the normally distributed random 

variable: coordinate of intersection of the weapon track 

extension with the vertical plane (a) through the 

submarine target center perpendicular to any plane con= 

taining both target and weapon at time t^ « (yd) 

DwTKt = minimum target depth to be considered (ft) 

MiN 

= initial acquisition range (yd) 

Q = number of depth subintervals in the depth interval 



Intelligence reports and estimates or CNO requirements may 

provide the target inputs. The weapon inputs are available from 

weapon designers at least in the form of educated estimates. The 

inputs RG and 0 are functions of the medium;> the weapon sonar and 

guidance systems, the target physical dimensions^^ and the orientation 

relationships between the target and weapon. In the absence of a 

specific theory these must be handled by simulation or sensitivity 

analysis. The input is a function of the capability of the 

weapon delivery system and the tactical situation as seen by the 

target commander. In the absence of specific considerations cri« 

terion one would determine to be zero. The input R. should be 

MIN A 


8o 


















selected close to the maximum weapon capability in accordance with 

criterion one. The function P may be very sensitive to changes in 

R and an examination of results over a range of values for R. may be 
A A 

appropriate. The input Q should be large enough to point up the 
effect of an irregularly shaped fuel expenditure rate curve but small 
enough not to excessively prolong the computer program. 
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CHAPTER V 


PROBABILITY OF KILL DAMAGE 

Consider a tactical situation in which a weapon explodes within 
a specified proximity of some portion of a submarine target and make 
the following assumptions* 

Assumption 1: The hit depth is a random variable which has an 
assumed probability density function* 

Assumption 2: The point on the target at which the hit occurs de» 
noted the hit point is a random variable which has a uniform distri® 
bution over the entire surface area of the target* The weapon 
approaches the target from astern, loses contact at some close-in 
range and continues on a straight track* Since the target may be 
in a turn, climb or dive, all points of the target are available 
as hit points. 

Assumption 3t The target operates in depth within specified minimum 
and maximum limits such that the maximum depth limit is not greater 
than the weapon maximum operating depth* Divide the target surface 
area into P area elements, each of which has a constant hull 
vulnerability throughout. Divide the interval of possible hit 

depths into Q equal depth subintervals such that the depth at the 

I 

lower boundary of the q th subinterval is 


\lN ^ 
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where 


^TM " target depth (ft) 


MIN 


= minimum target depth to be considered (ft) 


Define: ^ j^kill damage given an explosive weight w^J 

f„ (x) = probability density function of the random 
variable hit depth 

Pg^ = P ^hit depth is in the q th depth subinterval] | 
q = 1, Q 

Pg = P [^kit point is in the p th surface area element]| 


p 1, 


P 


HVp^(w 3 ) = hull vulnerability for kill damage of the p th 

target surface area element at the midpoint depth 
of the q th depth subinterval given an explosive 
weight w^. This number is one if kill damage occurs 
under the specified conditions and is zero if kill 
damage does not occuro 
W 3 = explosive weight (lb) 

III _i I I 

W , W = lower and upper limits respectively on the range 
of explosive weight (lb) 

Assuming as in Chapter IV 




^“im" “mih^ 
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Then 


P = ^ ^ a = 1 0 

% q 2 > ^ ^ 

p _ area of p th surface area element ! p s 1, »'», p 
Ep ~ total target surface area 


Pg = C''^) = 


HV (w.) P„ P„ 
pq ' 3 Ep Eq 


p=l q=l 


I V ff 19 


The table of P versus w- as varies from W to W 
is the input to the dynamic programming modelp 


Input Requirements 

1. Target 

Pgp = P l^hit point is in the p th surface area elementj 
p = 1, * • P 


= maximum target depth (ft) 

TM 

2. Weapon 

III _iIt 

W , W = lower and upper limits respectively on the range 
of explosive weight (lb) 

3o Target and Weapon 

HV (w^) = hull vulnerability for kill damage of the p th 

pq 3^ 

target area element at the midpoint of the q th 
depth subinterval given an explosive weight (w^) 
k. Miscellaneous 

D = minimum target depth to be considered (ft) 

MIN 

Intelligence or CNO requirements may provide D_„, P and some 

TM Ep 

of the information for HV (w-) . Weapon designers may estimate 

• I I j 

W and W and weapons engineers may provide extrapolated 

estimates of HV (w-). 

pq 3 ' 
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APPENDIX A 


FLOW CHART NOTATION 


The following notation is used in the flow charts of the remaining 
appendices. 


Function Box 


A — 


The function box contains a number or expression to be evaluated 
A which is to be placed in memory cell B* 

Choice Box 



The choice box contains a comparison of the number or express 
sion A to the number or expression B, Exits are chosen as follows® 

Exit is chosen if 
A> B 
A > B 
A = B 
A <B 
A < B 


Exit Label 


< 

< 
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In Box 


INPUT 


The in box contains an input to the program« In appendices 
C 1 through C 8 the output from one appendix is carried over intact 
into the next appendix. Inputs to be used for the first time in the 
sequence of flow charts are placed in a separate in box. 

Out Box 

OUTPUT 


The out box, one for each appendix, contains the output of the 
appendix program. 


Index Box 


1 i 


} 


The index box contains a manipulation on an integer in an 
index register and is the tool which provides for loops in the 
program. 
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Subroutine Box 



The subroutine box indicates a return jump to the subroutine as 
numbered by k. The flow chart of each subroutine is shown after 
the flow chart of the first routine using it* 

Return Box 



A return box indicates the end of a subroutine implying a jump 
back to the main routine. 

Fixed Connector 



A fixed connector is used to connect remote parts of the flow 
chart with one another. 


Assertion Box 


Note j 


I 

An assertion box refers the reader to a note found at the end of 


the flow chart. 
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Instruction flow is generally from left to right and from top to 
bottom. All boxes are entered from the left or from the top and all 
exits from boxes are made down or to the right. The arrow is used to 
further indicate the sequence of computer instructions. Arrows are 
used at all line intersections to indicate entry into the program 
between two boxes. 
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APPENDH B 

THE ALLOCATION MDDEL 
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Note 1: Provision is made for a set of payload weight limits (w ) in 

n 

case it is desired to examine a range of values. The input 

W„ is to be the greatest of the W , 

N ° n 

Note 2: The subroutine RJl computes the h^(w) and the Wj^(w) for 
i=l, 3* hg is handled separately. 
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Note 3 


Note k 

Note 5 
Note 6 

Note 7 
Note 8 

Note 9 


This loop provides that, for the W considered, the function 

is maximized over the . For hj^ the maximum meaningful value 

in both loops is subroutine is stopped by the 

comparison W ; D, However for h^, the maximum meaningful 

value in the two loops is different, the one in this loop being 

w^_ and the one in the other loop being W.,,„ s D. 

3b MAX 

This loop stores the h^(w) and w^(w) obtained as W is allowed 
to vary over its range of meaningful values« 

This loop maximizes the over j for a given k. 

The maximum meaningful value of w^ for a given W is 


min 


{”2K’ "} 


The h^Cw), WgCw), and W 2 *(w) are stored for each value of Wo 

This loop works backwards through the stages to retrieve 

the optimum weight values for each input payload maximumo 

The are the optimum volume sums to be compared with the 

input payload maximum allowed volume Vo If the solution is 

not optimal with respect to the restraints^, a new X is chosen 

and the program repeated. If the final solution is obtained^ 

then fuel weight is + I^^ weight for added attack speed is 

H and explosive weight is P . 
n n 
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APPENDIX Cl 
AVAILABLE ECHO EXCESS 



9k 



























































































































Note 1: This flow chart neglects the effects of ambient sea noise and 
of volume reverberation. Both are considered to have negli=> 
gible effect on the interference background masking level. 

If it is desired to consider the volume reverberationj, it is 
= WS - 2Hj^ + 10 log m^ log ( - | - ) 

+ 20 log (look) + RD; k 1, •«», K 

where m(R ) = volume reverberation masking level at range 100k yards (db) 

' V k 
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WS = on^axis source level (db at 1 yard from transducer) 

= one way transmission loss at range 100k yards (db) 

10 log m^ = volume scattering coefficient (db) 

(j^)j^ = volume reverberation index of the transducer (db) 
c s velocity of sound in the medium (yd/sec) 

AT = ping length (sec) 

RD = recognition differential of transducer (db) 


The index (J^)j^ t>e approximated by 


= 200 ko) 


= -10 log ^3 \ ^ P '^S )d_rL 


k = 1, 


K 


where m = search rate of turn (deg/sec) 

j = off-axis angle in the vertical plane (deg) 

P = off-axis angle in the azimuth plane at time of 
pulse transmission (deg) 

bg, b^, b|^ = normalized pattern functions of 

the transmitting (l^ 3) nnd receiving 
(2p k) patterns in the vertical (1^2) 
and azimuth (3^ ^) planes 

_f)_^ = solid angle subtended by the sound energy beamo 

This in effect assumes a medium of infinite extent 
in the direction of the sound beamo 
Note 2% M(SN) = SN + 10 log A f-DR -h RD 

where M(SN) = self-noise masking level (db) 

SN = equivalent isotropic self-noise level at search speed 
and depth (db) 
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Note 3t 
where 


Note hi 
where 

Note 5; 

Note 6; 
where 

Note 7t 
where 

Note 8; 


A f = receiving bandwidth (cycles/sec) 
DR = receiving directivity index (db) 
2H = WS + TS=M(SN) 


H s maximum possible on-^axis range for acquisition (yd)o 
This neglects noise levels except the self=noise level 
and neglects the weapon enabling level. 

TS = target strength (db) 

= 20 log look + 100 ak 

a = attenuation coefficient (db/yard) and it is assumed that 
the transmission anomaly is zero. 

This loop finds the upper limit (IOOK) on the acquisition 
ranges to be considered. 

M(TN)j^ = TN + 10 log 6 f=Hj^ + RDj k^l^ K 

M(tn)j^ = target noise masking level (db) 

TN = radiated noise level of the target (db at 1 yard from 
transducer) 


D. 


0, = arc sin 
k 


WS 


300 k 


= the angle above the horizontal at which an arc of radius 

100k yards and of origin at search depth will inter= 

sect the surface of the water 

D „ = weapon search depth (ft) 

Ww 

(j ), = surface reverberation index of the transducer (db) 

= angle of attack (pitch) of the weapon (deg) 

This is an approximation described in Chapter k of refer= 
ence k • 
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Note M(RS)j^ = surface reverberation masking level at a range of 

look yards (db) 

10 log m^ = surface reverberation scattering coefficient (db) 
Note 10: This subroutine adds the effects of two noise levels in 

decibels. A rounding off process is employed based on Figure 3 
of reference 3* A difference in levels of more than 10 db pro° 
vides a sum effect equal to the larger of the individual 
levels. Two equal levels have a sum effect equal to 3 
higher than either individual level. 

Note 11: ~ weapon enabling level (db). Values out to a range 

of lOOK* yards are inputs where is estimated such 
that it will be larger than the K found in the program. 
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APPENDIX C2 

SCANNED VOLDIffl RADIALS? GENERAL CASE 


INPUT 


INPUT 



PRRg f « » «• p D, 


TM 


APPENDIX Cl OUTPUT j- 


200 ^ 1 —^ i^ —^ (i®*!) 


P 


Note 1 




- ( 1 » 


r> 


1 - 1^1 






Note 2 



RJ3)-<l-^k 


P j2k * ^ 



^ 1 •» k-^kV 





^4k» “Jk‘ 

p*. 6, ® , D^, D^„, PRR^ 
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(RJ4] 


• ( 1 ^0 — { 


i 


Pjak • Pja(k+1) 


I^Note 91 


- { 1 » k ->-k ) — ( 1 » n -a»»a ) — 


P 


Note 10 


Pj»k-^Pj(k.n)H-l^k ; c/2Fm^ 



100k-^R 


J Ck»n) 


J(k=n) * "sln'(jV)c®® 



(D^ = Dv^)/3 


R. 


sln(j+|i)c©8p^^jg^^^ ^ j (k-a) 


p ♦ 

> 

®Us/3 

Jlk-n} 8ln(J+ji)co8pj^j^^^^ 

' ---- ^ 





(Note 12 


Pj(k.n) • P’ 


< 



(TTi 

:^k) 


k - n -^K, 



(RTN 


Pjl(k-n) • Pj2k 


< k » l^k^ - (n 2-^al - 


<D 



- (n - 1 -^n) - (2-^111 ) - 


Note 13 


‘- Mk - l^T) - (2_ ^m ) - (n » 2-^ny - 


<D 
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which meet the requirement that the echo will be received 
in the interval - «9 •> 


Note 2; The overlap in plane zero is computed first to provide the 
angle p* . Then p* = & and p’ define the limiting angles 

P jmk- 

Note 3; P ) for the i th P . 
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Note Uj 

Note 5: 

Note 6% 

Note 7: 

Note 8; 

Note 

Note lOj 
Note 11: 


After the computation of which indicates that the maximum 

acquisition range in plane j and all of the p , in plane j 

jniK 

have been computed^ the subroutine endso 

This indicates that there is no acquisition probability in the 
plane j and changes the range of planes to be considered from 
J to j\ \ 

In this case the p limited by the P ® and would have 

been larger if the overlap between pings had not been removed 
from the considered scanned areao 

This loop starts with the leftmost angle ( P® ) 

to be considered and computes each angle to be considered 
indexed by i and incremented by ^ to a maximum angle Pj^p® o 
This point is the commencement of the loop which will compute 
the P ,, and R., starting with j ^ =Jo 

JK JK 

There may be several equal P because the scan azimuth limits 
were reduced by the overlap <> The first of the equal p to be 

maximum ko 
yards/pingo 


considered will then be the one of maximum k 

C 


The ping interval range is 


2PRI^ 


The range at which the radial at angle Pintersects the 
surface of the water for j >> « p-is computed below* 

See Figure 20o The distance OC is desired* The search depth 
is (dWS/ 3) yards* 


length (CB) = / 3 sin(j^^p) 
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Figure 20. Computation of radial length for angle 

The plane j for which the scan outline is shown 

is perpendicular to the vertical xy plane. Angles 

CBO and DBO are right angles. Length A^ Ois %fee 

search depth of yards, 

l 
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Note 


Note 


length (OC) = length (0B)/cosp^j^ " D^^/3 sinCjHf-ji) eosp^j^ 


For j “ “JA the development is the same except that (j + }a ) 

is negative and the considered depth is the difference between 

the maximum target depth and the weapon search deptho 

TM 

12? If to the left of the radial of maximum range K. then m ~ 1 

J 

and the following angle will have a larger ko 
If to the right of the radial of maximum range then m ^ 2 
and the following angle will have a lesser ko There will be 
one or two radials of maximum rangeo That is 




13s If = Pj 2 K/ program 


program includes it, 


6 o Otherwise the 


lOU 







4PFENDIX C3 


SCANNED VOLUMESI GENERAL CASE 



Note 1| required iu the program but it may be desired 

to print it here as information. 
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APPENDIX C4 


SPECIAL VOLUME ELEMENTS* GENERAL CASE 
INPUT INPUT 


©-^ 


®TS» ^WD' 




W 


*Ot® It 


M 


APPENDIX C3 OUTPUT 


} 


^||360/t.|| 


n 


1 + 1 — 


1 ▼ X n 



°W5^0t^^ *^DS ^WD~^^OSD 


□ 


[(W'-w_„)/w. 


OSD''"It 


PRRgAj + l^U 


-lES 




T)- 




t V 


KB-^q> 



'^RJ6 


D« 



kB^q> 


OUTPUT 

LBDV. 


du'’-du''-du"” 

V^u? d = D+l, ***p DS* u = Ip 
BW^^l d = D+lp •“•p D"? u = Ip 

<S =9 

rr m ^ ^ ^ ood t\<\ c ^ o 


1« D^ u 


ASV, 



L 



= 1 » 

0 0 0 n? ® 
9^9 

U : 

^ 2, 

“°*P u 

qSV_ p* q : 

= Ap : 

B 9 u = 29 

0 0 

% u 


^Ad» 

®AId 

9 d s 

I9 

0 0 e 

9 

D» 

®A0d 

, tan^iod* 

d = Ip * 

e 0 

9 

D 


*^Bd» 

®Bid 

9 JfJ 9 ^ “ 

1 , 

000 

9 

DO 

®B0d 

• ^"*B0d5 

d - Ip ° 

0 0 

9 

D 

c=> 


Pjk’ 


i = 

J«p ““•p 

JO D 

® k s 

«n 0 0 e 

X 9 9 

^DS» 

®WD» 

®TS’ 

'^OSD® ^It* 

M 9 u 


P'» 


°WS» 


B 
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Note 10 
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["Note 13 


Rjl;^ 


A^tan(6/2) 



(RTN 



RJ14)_L 


I'Note 14 


(1/3) tan(6/2) C(D^ = E^)^P 



RTN 



I--- 


B - F^^IqDV 


du 



(RTN 


fRJ17>_L 


I" Note 11 


HL 


(TES ■= L^^^2)^(®'(<i=.l) “ 


(d-1) 


^d “ ^(d=l) 


HL. 


I— d 


— ( 1 ^r^r> — (HL^ = “ BS^)^CT-r 


^ ^(d-1) ’ ®^d^V 



(RTN 



depth. 

Note 2; The following comparison is used to determine the plane 
j to be used as a reference plane for the computations. 
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Note 


The following four function boxes set up reference constants 
for the subsequent computations« 

Note ^4-: The index d is used to identify the levelso The index (d^r) 
identifies the current plane j above the plane TJo The 
index (d-s) indentifies the current plane j below the plane 
T J- A T J o 

Note 5o This comparison determines if the level d is to be determined 
by a plane j> TJ or by a plane j<(TJ= ^TL)o 

It is only required for depths at or above the critical depths 
that is for d ^ D. 

Note 6: This subroutine computes the distance L^o See Figure 21* 

As shown ^ sin (TJ + 2 + p. ) > 

H(tj 2 ) (TJ-2 + ) and ATJ - 0. The distance to be 

computed is 

L, = length AB = length FG = length FE= length GE 
a 

Ld = Si 

Note 7: If the level length is longer than the distance TES 

traveled by the submarine during a complete weapon search 

turn, the critical depth is at a level above L, but below 
' a 

L/, ,N. Then the level defined by TES must be found and 
V““l J 

designated and the original level dropped from con^ 

sideration until later. If L is shorter than TES then it is 

above the critical depth and its parameters must be computed 

to find qSV , , the volume searched in the d th level on the u th 
odu' 

search turn, (q = A,B). 


112 











TJ 



Figure 21: A sample volume portion for computation of L.. The points 

d 

(Tji^ d) are the radial endpoints of the corresponding 


planes. The level is under consideration and 


Lj) is 

to 

be computed. 

The 

BC 3 - 

SC 3 

= length 

FE 

BS 3 - 

SB 3 

= length 

BG 

SS’ - 

SB 3 

= length 

HI 

SC^ - 

SC3 

= length 

HE 
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Note 8: See Figure 22. Subroutines RJ8, 9> 10 compute the distances 

and angles on the right portion of the figure which have not 
already been computed« Then subroutine RJll computes all of 
the special volumes in level L^o After RJll the level index 
as well as the index r or s^ is increased by lo If index s 
is changed it means that the plane which determined the level 
was above plane TJ« Thus d and s must be increased together 
so that the difference is constant« It is their difference 
that is used for the comparison described by note 5« 

Note 9r This subroutine computes the scanned volume in level 

(where d^D) on each search turn after the first <> The first 
comparison is used to find out if the u th search turn is the 
last effective search turn in level L^o If it is not^ then 
subroutines RJl2, 15/ ^nd 16 are used to compute the volumesc 
The cross-section in this case looks like that of the right 
portion of Figure 22. If the u th search turn is the last of 
the effective search turns in level L,^ then the normal cross^ 
section is cut at some point as in the left portion of Figure 22o 
In this case, the point where the cross-section is cut is 
found by one of the subsequent three comparisons and the cal^ 
culation of scanned volume is modified accordingly« 

Note 10: This corresponds to a vertical through the transducer which 

goes through a point in level such as point S in Figure 22o 
Note 11: This is used to compute the height of the critical depth 

above the search depth (length OR in Figure 22) and a correc¬ 
tion factor to be used to find the constants for level 
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m;. 








Figure 22: Vertical cross»section through levels and 

The right portion depicts the volume scanned on the first 
search turn. The vertical OAB extends down to the trans¬ 
ducer . 

Volumes corresponding to depicted areas are 
Area Volume 


HID 

DEIJ 

EJK 

CFG = EJK 
ABCF i DEIJ 
HDEK 
ABCG 

The distances used to 
length NQ = TES 
length BD = 
length BE = 
length AH = ®qid 
length AK = B^^^ 
length HK = 


LqDV 

MqDV 


dl 


dl 


RqDV 


dl 


RqDV^^ i RqDV^ 
MqDVdu i MqDV 


dl 


qSV 

qSV 


odl 


odu 

compute volumes in level L, are 

d 


length HI s b 


Id 


length JK s length FG - b 


od 


length AB s h^^ 
length AR = HL^ 
length BR = 
length DE = 
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Note 12; This is the volume computation for BqDV^^o This computation 
begins by turning the volume upside down so that level L. is 
above level Then the computation is the same as that 

for (corresponding to area ABDH in Figure 22) required 

later. 

Note 13; See Figure 23. 

Note ih; See Figure 2^. 

length A I = D tan0 
area UK = tan (6/2) 

volume AJKI = (l/3) (area IJK) (length Al) 

= ( 1 / 3 ) tan ( 6 / 2 ) tan 0 


Note 15; 


Note 16 ; 


Note 17: 


Note 18 ; 


The volume computations for levels below the critical 

depth are begun at this point. 

The volume computations for the last level use for the 

previous level and the angle j + d ■= r + p because there is 

no plane j on which to base the calculations. 

The calculations below search depth are Identical to those 

above after the H. below are substituted for the H, above. 

J J 

The final section computes the scanned volume on the u th 
search turn separately for above and below search depth. 
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Figure 23: Volume computation of a triangular 






EB 


= A, Length BC = A tan (h /2). Length ED = h^^ 
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The altitude lengths are AI and AH 













APPENDIX C5 


FIRST TURN SCANNED VOLUME: i TH PING 



Note It p is now the angle between and the radial 
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APPENDIX C6 


SCANNED VOLUMES SUBSEQUENT TO FIRST TURN: i TH PING 
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0 — 


CH 


Note 7 


> 

{Note 8 


1-- 

1 


<D 


(CH)2tan(&/2) - ^®ql(d-l) " = CH)tan0qj^ 

- (l/3)tan(6/2)tan^qjj [(CH)^ - (B^j^ = (u-DTES)^ 


+ ASV^ ASV 



^ ^ ) -[£^Td ^- 


Note 9 




1 < f Note 10 

1 ^ 

' J_1_ 

~n[> 

L|B^^|0Tu-2)roS~TcH] 

^ fNote 11 



AST ^ qSTpj^^ZsT]- 0 


(l/3)tan(&/2)tan^qjj [(CH)^ - = (u=2)TES){| 

= 9M7^u * ^SV^ASV 






^Note 11 


l> 



n 


(CH)^tan(&/2)h ^ + aSV-^ASV 


IB _ - (u-l)TES : CH 


,< 

'Note 12 

V 

1 



•0 


(CH)2t<in(b/2)[h^^ - (B,i(a.i) - (u-l)TES - 

- (l/3)t«n(&/2)tan^'^ jj [(CH)^ - (B^jj - (a-l)TES)^] ♦ AST-5-AST 


©-C 


1 + d ->d 


1- asv ->■ - (1 » u -$> u 




. f 


Note 14 


I"Note 13 


0 


q : A>=-L 


-® 


- <B-^q 


D ->d. 


D«->d, 


0 
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Note 1 
Note 2 

Note 3 


Note 4 
Note 5 


Note 6 

Note 7 

Note 8 


P is now the angle between and the bisector of the ith pingo 
: In this case the target volume cylinder is large enough such 

that the scanned volumes after the first search turn are not 
effected, 

: In this case the target volume cylinder is small enough such 

that some of the scanned volumes after the first search turn 

must be modified,, This is accomplished above and below the 
search depth and for each search turn by summing the modified or 
unmodified volumes over the levels to obtain First 

each search turn above search depth is considered^ then each 
search turn below• 

In this case the original is not changed« 

: In this case the original LqDV^^ and MqDV^^ are not changed 

but RqDV^^ is modified and must be recalculated,, The 
original cross-section appears in Figure 25• Using formulas 
of subroutines RJ13 and RJl4 the new volume P is P « volume 
FGED + volume ADEC - volume ABGF« 

: In this case the original is zero^ the original LqDV^^ 

is unchanged and the modified MqDV^^ is computed from the 
formula for subroutine RJl3« 

: In this case the original MqDV^^ and are zero and the 

modified is computed from the formulas for subroutines 

RJ13 and l4 as the difference of two volumes* 

In this case the contribution of level L, on the u th search 

cl 

turn is zero* 
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Figure 25: Computation of modified volume RqDV^^o Length FG is 

®qO(d-l) ' Area BGH represents the original 

volume RqDVj^. Line AF is part of the vertical through 
the weapon. 


125 

























Note 

Note 

Note 

Note 

Note 

Note 

Note 

Note 


9: The contributions of levels above the critical depth have 

now been determined • 

10s The contribution of level is not changedo 

11: The missed volume on the previous search turn has been reduced 
leaving less to be picked up on the unchanged second turuc 
The formulas are those of subroutines RJ13 and RJl4^ 

12: The missed volume in both the current and the previous 

search has been reduced» The formula is identical to that 
discussed in Note 7« 

13: The scanned volumes indexed by one set of u^ and i are 
now computedo The next step is to change u and compute the 
next set of volumes• 

14: Scanned volumes on the i th ping of all search turns above 
search depth are now computed« The next step is to compute 
the volumes below search depth« 

15: The increments of scanned volume for ping i^ 
computed. The next step is to compute 

16: Having found the index is changed to find 

(i + l)u 
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APPENDIX C7 


TARGET volume: i TH PING 
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Note li The problem is to compute areas as depicted horizontally and 
vertically lined in Figure l6o Using Figure l6^ with origin 
at point C and x axis toward point the equation of the 
target volume cylinder circumference on ping i is 
(x + R )^ + = R.^ 

Changing to polar coordinates (r^ p) with the same origin 
and measuring angles positive counterclockwise from the 
X axis^ the equation becomes 

+ (2R cosB ) • r=(R? =R^) - o 

Solving for r, 

r = R^ ^ I - (R /r. sin =R cosP 
i 'o'l' o 

To find the area for ping i, r(r d p ) is integrated from 
p' + (i - 2)6 to p*+ (i - 1)6 . 


Thus 

area = 

Let p' + (1-1)6 = B 
p» + (1-2)6 = A 


p'+(1-1)6 . 

dp 

p'+(1-2)6 


Then 
■B 
A 


r r2dp= r [(Ri^. 2) ^ 2 2^^^2p 

- 2 RoR^Vl “ (R^/llj^) Sln^p cosp 

J (R^^ - R^^) dp = (R^2 - R^^)(B - A) 


dp 
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2 R ^ cos^B dB = 

O o 


p + 


(1/2) Sin 2 p 


B 

A 


Let X = sinp 

Then dx = cosp dp 


- 2 R^R^ J* - (RqA^)^ sin^p cosp dp 

' ■ ^ * 0*1 fk 


2^2 dx 


Let sinp = (R^/R^) x 
Then cosp dp = dx 


" ^ / \/iT7r~/R^ 1^dx = - 2 Rj^ 


2 ^ 

COS p d p 


= » R^^ jp ^ (l/2)sin2p 


Therefore I 



r^ dp = (R^^ - R^^) (sin2B - sin2A) / 2 
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APPENDn C8 


ACQUISITION PROBABILITIES! I PINGS 


. INPUT _^ , INHJT-- - 

1 ) - ( APPENDIX C7 OUTPUT ) — ( Av, W% W» ) - (1— 



w 




SVi„ P«lATi-^8i„ 

—(1 + i-^i )—(i : 


1 u-^u 





' -( 1 + ss ) —(1—>-u)- 

0 —>g»j„T 

—^<l->i> 

\ _/ \ -J 

®MU 

\ ... / 


^gju^gsl ^ ^ •*• — ( 1 ^• i-^i' 


L ( 1 ♦ g -^u 




® 


(w*/Aw) “ ^’/aw) —^ R 



1 


W* —> w. 




[(W^ - W^gp) PRRg/W^^ 


-^l 


l^AP, 


cSsS 




> - 

- --- 

-(l^i>JL 

(l-gljAP^-^APj^ 

—( 1 + i-^i \ —Ti : i}^ 




-{1 * r-^T~)- 


Aw + w. 


l(r=l)‘^''lr 




<D 


OUTPUT 


“ir' " 
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APPENDn D 


PROBABILITY OF HIT 
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Note Ij The problem constants are now determined. Next the 


probability P„.. is determined for 

HJK 


j = l, J,jk = l^ 


K. 
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